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(57) ABSTRACT 

Shape of a multi-core optical fiber is determined by position- 
ing the fiber in an arbitrary initial shape and measuring strain 
over the fiber’s length using strain sensors. A three-coordi- 
nate p-vector is defined for each core as a function of the 
distance of the corresponding cores from a center point of the 
fiber and a bending angle of the cores. The method includes 
calculating, via a controller, an applied strain value of the 
fiber using the p-vector and the measured strain for each core, 
and calculating strain due to bending as a function of the 
measured and the applied strain values. Additionally, an 
apparent local curvature vector is defined for each core as a 
function of the calculated strain due to bending. Curvature 
and bend direction are calculated using the apparent local 
curvature vector, and fiber shape is determined via the con- 
troller using the calculated curvature and bend direction. 

20 Claims, 2 Drawing Sheets 
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SHAPE SENSING USING A MULTI-CORE 
OPTICAL FIBER HAVING AN ARBITRARY 
INITIAL SHAPE IN THE PRESENCE OF 
EXTRINSIC FORCES 

5 

STATEMENT REGARDING FEDERALLY 
SPONSORED RESEARCH OR DEVELOPMENT 

The invention described herein was made by employees of 
the United States Government and may be manufacmred and 1 0 
used by or for the Government of the United States of 
America for governmental purposes without the payment of 
any royalties thereon or therefor. 

ORIGIN OF THE INVENTION 15 

The invention described herein was made by employees of 
the United States Government and may be manufacmred and 
used by or for the Government for governmental purposes 
without payment of any royalties thereon or therefor. 20 

TECHNICAL FIELD 

The present disclosure relates generally to fiber optic shape 
sensing, and in particular to a system and method for fiber 25 
optic shape sensing using a multi-core fiber having an arbi- 
trary initial shape and extrinsic forces acting on the fiber. 

BACKGROUND OF THE INVENTION 

30 

Single-core and multi-core optical fibers are used in a 
variety of applications. In a single-core optical fiber, a single 
light-guiding core is contained within reflective cladding 
material. A multi-core optical fiber includes a plurality of 
such cores. Strain sensors such as Fiber Bragg Gratings 35 
(FBGs) may be attached to or embedded within the fiber 
along the fiber’s length. FBGs may be formed by laser-in- 
scribing, writing, or otherwise embedding a periodic varia- 
tion of refractive index into the cores of the optical fiber. This 
effectively creates an in-line optical filter designed to block 40 
particular wavelengths of light transmitted through or along 
the core. Alternatively, Rayleigh scatter detectors can be used 
to detect elastic light scatter occurring within a given core at 
specific axial locations of the optical fiber. Using these and/or 
other strain sensors, the bending geometry of the optical fiber 45 
can be calculated. The calculated geometric data may be used 
in various ways to approximate the shape of an object to 
which the cable is attached. 

Strain information from axially co-located strain sensors 
can be used to estimate the bending at each co-located sensor 50 
location along the fiber. Common methods for accomplishing 
this exist, including, those outlined by S. Klute et al. in “Fiber 
Optic Shape Sensing & Distributed Strain Measurements on 
a Morphing Chevron”, 44 ,,! American Institute of Aeronau- 
tics and Astronautics (AIAA) Aerospace Sciences Meeting& 55 
Exhibit, #AIAA 2006-624, Jan. 9-12, 2006, Reno, Nev. Such 
methods are highly dependent on the accuracy of successive 
strain measurements used to calculate the cable’s bending 
parameters at discrete segments of the cable. 

When a multi-core cable in particular is subjected to bend- 60 
ing, the strain in each core depends on the curvature and 
direction of the bend, as well as on the arrangement of the 
various constituent cores within the cable in relation to the 
bending direction. In accordance with linear elastic tube 
theory, a light-guiding core positioned on the inside of a bend 65 
experiences a stress, i.e., a negative strain, while a core posi- 
tioned on the outside of the bend experiences a positive strain. 
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The amount of strain is proportional to the bend radius and to 
the position of each core relative to the center of the bend 
curve. Therefore, multi-core fibers can have additional utility 
in comparison to single core fibers when used for structural 
shape sensing and end effecter tracking. 

Conventional fiber optic shape sensing approaches initial- 
ize and calibrate the sensors while the fiber is straight, i.e., 
when the fiber initially has zero curvature and zero strain. 
U.S. Pat. No. 7,813,599 to Moore (hereinafter Moore ’599), 
which is hereby incorporated by reference in its entirety, is 
one such approach. However, in practical applications it may 
be relatively difficult to maintain a known core orientation, 
and thus the angular core orientation is not always known. 
Using conventional methods, the calibration of strain sensors 
while the fiber is in an arbitrary initial shape, which may or 
may not be straight, may introduce a bias error that can lead to 
a substantial error in shape measurement. 

Additionally, conventional shape sensing approaches do 
not adequately account for the effects of common extrinsic 
forces on the fiber such as twisting and/or stretching. Extrin- 
sic axial loading may be introduced into the fiber through 
attachment mechanisms such as adhesives or tape, by general 
handling of the fiber as in a catheter application, or by friction 
inside of a sheathing winch may house the fiber in a given 
application. Such extrinsic loading can cause shape measure- 
ment errors to occur in prior art approaches for determining 
the shape and/or position of a multi-core fiber. 

SUMMARY OF THE INVENTION 

A system and method are disclosed herein for determining 
the shape, and thus the end position, of a multi-fiber cable 
having an arbitrary initial shape, i.e., a shape that is not 
necessarily straight. For instance, the cable may be curved, or 
the cable may be wound on a spool in a helical shape. The 
fiber of the cable may be free-sleeved as in Moore ’599 or it 
may be bound to a sleeve of the cable without changing the 
present approach. Moreover, the present approach may be 
used even when extemal/extrinsic pinching and/or pulling 
forces are present, which may ultimately cause extrinsic 
twisting and/or elongation of the fiber. As noted above, such 
extrinsic forces are not adequately accounted for in the prior 
art. 

Additionally, the present approach enables the accurate 
calculation of curvature and bend direction using a multi-core 
fiber with an arbitrary angular core orientation. Prior art 
approaches require a carefully calibrated angular core orien- 
tation before any shape measurements can be made. The 
present approach by contrast enables initialization of a multi - 
core fiber with arbitrary core orientations in an arbitrary 
shape. This in turn enables the use of substantially longer 
lengths of fiber in practical applications. In the absence of the 
present system and method, initialization of a long length of 
fiber would require a complex apparatus and/or a large work 
area. In addition, the present curvature calculation approach 
is insensitive to external strain-inducing effects such as tem- 
perature and tensile stress. 

The present method also allows for shape sensing of com- 
plex structures without calibrating angular core orientations 
relative to the structure. As an example if a sleeve is woven 
into the fabric of an inflatable structure for the purpose of 
containing a multi-core shape sensing fiber, it can prove dif- 
ficult to maintain a known angular core orientation while 
insetting the multi-core fiber into the sleeve. The present 
approach eliminates the need for angular core orientation 
cognizance during or subsequent to such an installation. 
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In particular, a method of determining the shape of a multi- 
core fiber optic cable includes positioning the cable in au 
arbitrary initial shape and measuring strain over a length of 
the cable using a plurality of strain sensors while the cable is 
held in the arbitrary initial shape. The strain sensors are dis- 5 
tributed along a corresponding length of each of the cores of 
the cable. The method may include defining a three-coordi- 
nate p-vector for each core as a function of the distance of the 
corresponding cores from a center point of the fiber and a 
bending angle of the corresponding cores with respect to a 10 
local axis of the fiber, wherein the coordinates of the p-vector 
define a plane. 

The method may also include calculating, via a controller, 
an applied strain value of the fiber using the p-vector and the ; 5 
measured strain for each core, and calculating strain due to 
bending as a function of the measured strain and the applied 
strain value. Additionally, the method may include defining 
an apparent local curvature vector for each core as a function 
of the calculated strain due to bending, calculating curvature 20 
and bend direction of the fiber using the apparent local cur- 
vature vector and an initial curvature vector, and determining 
the shape of the cable via the controller using the calculated 
curvature and bend direction. 

A system for determining a three-dimensional (3D) shape 25 
of a multi-core optical fiber includes a fiber optic cable and a 
controller. The cable includes an optical fiber with light- 
guiding cores and strain sensors. The strain sensors each 
measure strain at a different axial location of a particular core. 
The controller calculates, using the measured strain, an 30 
applied strain value as an axial strain component experienced 
by all of the cores due to at least temperature and axial load 
acting on the cable. The controller also calculates strain due to 
bending of the fiber as a function of the measured strain and 
the applied strain value, and defines an apparent local curva- 35 
ture vector for each core as a function of the calculated strain 
due to bending. 

Additionally, the controller calculates curvature and bend 
direction of the fiber using the apparent local curvature vector 
and the initial curvature vector, generates a curvature func- 40 
tion, a bending function, a twisting function, and a strain 
function by performing a curve-fitting, operation, and calcu- 
lates a function of torsion vs. fiber length of the cable, includ- 
ing differentiating the bending function, and then adding the 
twisting function. The controller can determine the shape of 45 
the cable in part by solving a set of Frenet-Serret equations 
using each of the curvature function, the bending function, the 
twisting function, and the strain function. In this manner the 
controller can derive the shape of the cable in a three-dimen- 
sional space. 50 

In another embodiment, a method of determining the shape 
of a multi-core fiber optic cable in the presence of extrinsic 
twisting and stretching forces includes recording discrete 
measured sets for each of: curvature, bending, twisting, and 
strain of the cable, and then performing a curve-fitting opera- 55 
tion on each measured set to generate a curvature function, a 
bending function, a twisting function, and a strain function. A 
controller is then used to calculate a function of torsion vs. 
fiber length of the cable, including differentiating the bending 
function and then adding the twist function. The controller 60 
determines the shape of the cable in part by solving a set of 
Frenet-Serret equations. The equations use each of the cur- 
vature function, a bending function, a twisting function, and a 
strain function to thereby derive the shape of the cable in a 
three-dimensional space. 65 

The above features and advantages and other features and 
advantages of the present invention are readily apparent from 


4 

the following detailed description of the best modes for car- 
rying out the invention when taken in connection with the 
accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a schematic perspective side view of a system for 
determining the three-dimensional (3D) shape and end posi- 
tion of an optical fiber. 

FIG. 1 A is a schematic cross-sectional illustration of a fiber 
optic cable used in the system shown in FIG. 1. 

FIG. 2 is a flow chart describing an example method for 
determining the 3D shape of an optical fiber using the system 
shown in FIG. 1. 

FIG. 3 is a schematic illustration of a unit circle for defining 
local axes and core locations. 

FIG. 4A is a schematic illustration describing the specific 
case of a fiber having a symmetrical three-core geometry. 

FIG. 4B is another schematic illustration of a fiber having 
a symmetrical three-core geometry. 

FIG. 5 is a schematic illustration describing the general 
case of an example M-core fiber cross-section. 

FIG. 6 is a flow chart describing an example method for 
determining 3D shape of the example cable of FIGS. 1 and 1A 
in the presence of bending, stretching, and twisting forces. 

DETAILED DESCRIPTION OF THE INVENTION 

Referring to the drawings, wherein like reference numbers 
represent like components throughout the several Figures, an 
overview of the basic construction of a multi-core fiber optic 
cable 1 0 is presented first with respect to FIGS. 1 and 1 A. This 
description is followed by an explanation of how a controller 
26 may be used to execute the present method 100, an 
example of which is shown in FIG. 2, in order to calculate 
strain and bending data for the cable 10 when initialized from 
an arbitrary initial shape, which may be straight or non- 
straight. Additionally, the controller 26 does not require an 
initial zero-strain state to function, and thus overcomes limi- 
tations of prior art systems in this respect. The handling of 
extemal/extrinsic forces on the cable 10 is discussed below 
with reference to FIGS. 5 and 6. 

The example cable 10 of FIG. 1 has a longitudinal axis or 
centerline 11 and a length (L). A plurality of light-guiding 
cores 12 are arranged within the cable 10. A three-core 
embodiment will be used hereinafter for illustrative consis- 
tency, although the present invention is not limited to this 
particular embodiment. Bach core 12 is positioned with 
respect to the centerline 11 and is surrounded by a layer of 
reflective cladding material 14. The cores 12 and the cladding 
material 14 may be constructed of a suitable polymer. The 
cores 12 have a high transparency, while the cladding material 
14 has a low refractive index relative to the material of the 
cores 12. Thus, any light emitted into one end of each core 12 
is efficiently transmitted along the entire length (L) of the 
cable 10 via total internal reflection, as is well understood in 
the art. 

Structural shape sensing and end-effector tracking may be 
performed using multi-core optical fibers as noted above. 
Strain information from axially co-located sensors 24, such as 
Fiber Bragg Gratings (FBGs) or Rayleigh scatter detectors, 
can be used to estimate the bending at each co-located sensor 
location along the fiber. Because the strain measurements are 
referenced to some initial state, prior work such as that of 
Moore ’599 can be used to calibrate the sensors 24 while the 
fiber is in a straight state, i.e., a state having zero curvature. 
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Certain mathematical solutions presented in Moore may be 
used in the present approach, particularly those solutions 
using known Frenet-Serret curvature equations. Flowever, 
critical adjustments are made via the present approach for 
calibrating the sensors 24 in fibers having an arbitrary initial 5 
state, i.e., a shape that may or may not be straight. The present 
approach is therefore directed toward the calibration of the 
sensors 24 and a determination of the true orientation of the 
cores 12 while the cable 10 is positioned in an arbitrary initial 
shape. to 

Multi-core fiber optic shape sensing requires that strain 
measurements be taken in each core 12 in order to properly 
determine curvature and bend direction. Because such strain 
measurements are relative to some initial material state, the 
initial or reference state must be accurately defined. The 15 
initial state can be defined by the curvature and bend direction 
of the cable 10 for all measurement locations. When the initial 
state is straight, subsequent bent shapes will result in com- 
pressive strain for cores 12 on the inside of the bend and 
tensile strain on the outside of the bend. As the current state of 20 
the art relies on the assumption that the initial shape is always 
straight, any deviation from that assumption would tend to 
cause error in the measured shape. The present approach as 
set forth below with reference to FIGS. 2-6 helps solve this 
particular problem. 25 

Referring briefly to FIG. 1A, the cable 10 shown in FIG. 1 
includes an optical fiber 15 having multiple cores 12, the 
cladding material 14 noted above, and a protective polymeric 
layer or coating 21 bonded to the cladding material 14. The 
cable 10 may also include a protective sleeve 23 constructed 30 
of a suitably rugged polymeric material, with the sleeve 23 
circumscribing the optical fiber 15 along to form an interface 
25. The optical fiber 15 in one embodiment is unbound, i.e., 
free to rotate or twist within the sleeve 23 at or along the 
interface 25 with little to no friction, as indicated by the 35 
double arrow 4. That is, the sleeve 23 is not bonded, adhered, 
or otherwise attached to the optical fiber 1 5 at the interface 25 
anywhere along the length L, thereby allowing the optical 
fiber 15 to freely twist and untwist with respect to the sleeve 

23 as needed. In another embodiment the sleeve 23 may be at 40 
least partially bound at the interface 25, such as by bonding or 
attachment, or may otherwise experience surface friction 
between the sleeve 23 and the fiber 15. The present approach 
works well with both embodiments. 

Referring again to FIG. 1, the cores 12 may be positioned 45 
substantially equidistant from and symmetrically with 
respect to the centerline 11, although other configurations can 
be used within the scope of the invention. Each core 12 is 
concentrically surrounded by the cladding material 14. For 
clarity, the cladding material 14 is shown only at the end 50 
portions of FIG. 1 in order to show the cores 12 and other 
internal detail of the optical fiber 15. Likewise, the sleeve 23 
of FIG. 1A is omitted from FIG. 1 for further clarity in 
depicting the cores 12. 

The cable 10 of FIG. 1 includes a plurality of sensors 24 55 
that are inscribed, written, embedded, or otherwise provided 
at intermittent axial positions along the length L. The sensors 

24 are configured to measure a set of strain data, and to relay 
the set of strain data, collectively represented in FIG. 1 as Sjsb 

to a sensor interrogation device or controller 26. The control- 60 
ler 26 is configured to determine the shape of the optical fiber 
15 via the methods 100 and 200 as set forth below with 
reference to respective FIGS. 2 and 6 using the strain data S^. 

In an example of a potential utility for the cable 10, 3D 
shape and/or end position sensing of the cable 10 can be 65 
achieved in conjunction with a device 29, such as but not 
limited to a piece of minimally invasive surgical equipment or 
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a tethered device such as a submersible vessel used in deep 
ocean exploration, a tethered camera of the type used during 
on-orbit heat shield tile inspection or other types, tethered 
robots used for search, rescue, and recovery operations, e.g., 
in collapsed mineshafts or in structures having compromised 
stability after an earthquake, etc. 

The device 29 may be connected to an end 27 of the cable 
10, with the sleeve 23 (see FIG. 1A) of the cable 10 in turn 
connected to an object 39 such as another set of motion- 
control cables or wiring used for controlling the direction, 
operation, and/or motion of the device 29 as needed. The 
controller 26 determines the 3D shape and end position of the 
optical fiber 15 of FIG. 1A as set forth below, with the 3D 
shape and end position usable thereafter for a variety of 
purposes, including for determining the precise position in 
3D space of the device 29. 

In addition to the end 27, the cable 10 has a bound end 20 
providing a calibrated end position. In other words, the cable 
10 is secured at the bound end 20, such as by attachment to a 
stationary member 22, e.g., a surface of a surface, nautical, 
submersible, or orbital vessel, vehicle, or other device used in 
conjunction with the device 27. Flowever configured, the 
bound end 20 represents a known initial condition usable by 
the method 1 00 for determining a shape of the optical fiber 15, 
and maintains a known initial tangent vector and a calibrated 
twist orientation as explained below. The shape the cable 10 
takes with respect to the bound end 20 may be non-straight. 
For instance, a helical shape is common on a spool-wound 
cable. 

As noted above with reference to FIG. 1 A, the optical fiber 
15 may be unbound or freely-sleeved, i.e., allowed to freely 
twist within the sleeve 23 at, along, or with respect to the 
interface 25 as its own elasticity governs in order to minimize 
internal stresses along its length L, while also conforming to 
the turns and bends of the cable 1 0. Allowing the optical fiber 
15 to freely twist and untwist in this manner enables the laws 
of elastic tube theory to apply, while also permitting explic- 
itly-defined functions of curvature, bending direction, and 
torsion to follow the known Frenet-Serret curvature equa- 
tions. The application of such equations is set forth below and 
in detail in Moore ’599. However, the present approach is not 
limited to freely-sleeved fibers. 

Referring to FIG. 2, an example embodiment of the present 
method 100 is described for sensing the shape of an object 
from an arbitrary initial state or shape. 

Initialization 

At step 102, the cable 10 of FIG. 1 is placed into an 
arbitrary initial shape, which includes any shape that the user 
desires, whether straight or non-straight. For instance, the 
cable 10 may be curved, straight, or shaped as a helix. The 
cable 10 thus has a known or determinable initial curvature 
vector when positioned in the arbitrary initial shape, which 
may be recorded by or otherwise made available to the con- 
troller 26. The ability to use an arbitrary initial shape enables 
initialization of a relatively long length of fiber while the 
cable 10 is, for instance, compactly held on a spool. The 
actual shape of the cable 10 is not limited. 

As part of step 102, the fiber 15 of FIG. 1A may be con- 
nected to a suitable interrogating system such as the controller 
26 of FIG. 1 once the fiber 15 is in the arbitrary initial shape. 
The controller 26 is then initialized, a step which may include 
initializing strain measurement via the sensor 24. 

Strain Measurement 

At step 104, the fiber 15 may be placed in a shape to be 
measured once the controller 26 has been properly initialized 
at step 102. The controller 26 may then request measurement 
of strain via, the sensors 24. There can be three or more 
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sensors 24 at each sensor location along the length (L) of the 
cable 10 shown in FIG. 1, typically one sensor 24 per core 12 
per location. For instance, a controller 26 that is interrogating 
sensors 24 in the form of FBGs may record the Bragg wave- 
lengthforall sensors 24 in the cable 10. Othertypes of sensors 
may be used in lieu of FBGs to measure the strain in the cable 
10 . 

Referring briefly to FIG. 3, certain terms used in the equa- 
tions described below are shown graphically for illustrative 
clarity, i.e., the local axes x, y, and z the locations of a core 12, 
and the bend direction 0 4 . The measured strain is incurred by 
a number of factors. The fiber 15 of FIG. 1A may experience 
strain that is not due to bending, e.g., strain due to temperature 
changes along its length, tensile loading, humidity variations, 
etc., in addition to the strain that is experienced due to bend- 
ing. The additional strain not due to bending is referred to 
hereinafter as the applied strain. Due to linear superposition, 
the measured strain at a given core 12 can be determined as 
follows: 

e,=-K d t cosfe^-0,) 

where: 

€ r i is total strain measured at the i-th core; 
e a is the applied strain experienced by all cores along their 
neutral axes due to temperature, axial load, etc.; 
e, is axial strain at the i-th core due to bending; 
k is curvature (a scalar); 

d, is distance from the i-th core to center of fiber; 

0 4 is the angle from local y-axis to bend axis; and 
0 , is the angle from local y-axis to the i-th core. 
Calculating Core Strain Due to Bending 

At step 106, assuming the applied strain components are 
experienced equally by each core 12 of FIG. 1, the bending 
strain components can be calculated as follows. First, one 
may define a “p-vector” for each core 12 as a function of the 
distance d ( and angle 0 ,.: 


the center point of fiber 15 (see FIG. 1 A), and has magnitude 
equal to: 


5 




Sj 

dj 


Using unit vectors j and k in the local y- and z-axes, respec- 

io tivel y ; 


Kappj = - -2- (cos#,- j + sin #,-£) 


The vector sum of the apparent curvature vectors for M cores 
is then formulated as: 


20 


= ~Z f cos9 <7'-£ fsinft} 


( 1 ) 


25 The magnitude of this vector from equation (1) is: 


\ K app\ : 


30 


M M 

Zl cos9i + z 


m 

— sin#,- 
dj 


( 2 ) 


By substituting the expression for strain from equation (1): 


35 

M M 

K app = cos (#*, - #,-)cos#,-} + cos (9t, - #,-)sin#,-£ 

i-i i - 1 


' djcos 9j 
pi = djsin9j 
£ t,i 


The magnitude of this vector is calculated as: 


These vectors define the three dimensional (3D) coordinates 
for each core 12 where the first and second components are 
the y and z components of the core location and the third 
component is the measured strain in the core 12. The three 
coordinates define a plane. The nonnal vector of this plane is: 


The applied strain value (e a ) at y,z=0 is the axial strain 
component experienced by all cores 12 due to temperature, 
axial load, etc., and is equal to: 


n-p 1 n-p 2 n-p 3 (0) 

e “~ n-[ 0 0 1 ] “ n ■ [ 0 0 1 ] “ n ■ [ 0 0 1] 

The strain due to bending (e,.) can then be calculated as the 
total measured strain less the applied strain, 

e i =e ti -e a 


45 


\ K app\ : 


(M 


cos(9b - 9j)cos9j + Z cos (9t - #,)sin#,- 


M 

Zc 

Vi-i 


(3) 


In general, the core locations may not be symmetric about 
the fiber axis, so equation (3) can be rearranged to express k 
50 explicitly. Incorporating equation (2) for the magnitude of the 
apparent curvature vector gives: 


F'wl (4) 


/ ( M 

\ 2 , 

' M 



J 1 z cos(fy + 

>; < 

30S (9^ - 9,)sm9j 

\ Vi=l 

) \ 

<i+i 


) 


M > 

2 

M \ 

2 

A 

X I cos9i 

, i=l ; 

+ 

Z u in " 

V 1 = 1 


i( M 


i 2 

( M 

r 2 

J |^Z cos {9b - 9j)cos9j 

h 

|^Z cos (#£ - #,)sin#,J 


The strain due to bending is then used in subsequent steps. 
Calculating Curvature and Bend Direction 

At step 108, an apparent local curvature vector is defined 
for each core 12 , which points to the respective core 12 from 


65 In the case of non-symmetric cores 12, the bend angle must 
first be determined using any previously reported suitable 
method. For s ymm etric cores 12, further manipulation yields 
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a simplified result that does not require the calculation of bend 
angle. 

Equation (3) may be expanded as follows: 


\^app — 

17m 7 7w 

k . / 1 y (casOpCOsO, 4- sirA,sin#,)cos#, 1 + I £ (cos#,, cos#, + sin#,,sin#,) 

Distributing: 


10 

Curvature can be calculated using the measured apparent 
curvature vector magnitude from equation (2) above: 


5 _ 2\K app \ (6) 

K ~ M 

Bend direction is calculated by determining the angle of 
the apparent curvature vector relative to the local y-axis (in 
10 thiscase, a bend direction of zero means the curvature is in the 
direction of the y-axis). 


l^appl — K 


/ M M \ 2 

| cos#s Z cos 2 #, + sin#,,X sin#,cos#,J - 

( M M \ 2 

| cos#,, X cos#, sin#,- 4 sin#,,X sin #,- j 


15 


Using trigonometric identities: 


20 


tappl — K 


{ cos 9 b M sin 9 b » r 

l-j^ i 2 (l+cos2fll) + -^2 sin2#,j + 

tcos 9 b M sin#* » t 2 

I ■ sin2# ; + -jZ-Z (1 -cos20 ; )l 


25 


Simplifying: 


30 


\ K app\ — K 




+ f, cos20,-j + 2^ 2 sin2#, 
\m- f cos29,j) 


COS 9 b M . ^ 

— - 1 Sin2 e, 

2 i = 1 


sin &b( 

' — 


35 


When the cores 12 of FIG. 1A are symmetrically spaced 40 
about center of the fiber 15: 


2tt 


By identity, when equation (5) is satisfied: 


M 


^ COS20; 


M 


V 1 sin20; 


s K M ^) cos ( 20i+(M_1) ^) 


2 n 


sinj 

( 2n 
M — 
l M 

sinj 

( 2 7T\ 

[Mj-MM-U-J 

2tt 


Therefore: 



50 


55 


60 




M 2 , . . 

= K^j (cos 2 % + siirOb) 


kM 

~T 


65 


6 b = cos 


■ = sin j j- = tan 


(7) 


Correcting the Measured Curvature Vector 

At step 110, the measured curvature vector is determined. 
The measured curvature vector is based on strain values that 
are referenced to the initial fiber shape, i.e., the arbitrary 
initial shape set at step 102. Therefore, the curvature vector 
that is calculated in step 106 will not result in an accurate 
shape calculation. In order to correct the measured curvature 
vector, the curvature vector for the initial shape is added to the 
measured curvature vector, resulting in the corrected actual 
curvature K acl as shown as: 

K„P,= K„ pp +K lnUa i 

where K. tntadt depends on the arbitrary initial shape, e.g., zero 
for a straight shape or some non-zero value for a helix, curve, 
or other initial shape. The actual curvature vector can then be 
used to determine curvature and bend direction as is shown in 
equation (6) and (7), respectively. In the non-symmetrical 
case, equation (6) is replaced with equation (4). 

Calculating Shape 

At step 112, the shape of the fiber 15 of FIG. 1A to be 
measured is then calculated by the controller 26 of FIG. 1 
based on the corrected curvature and bend direction measure- 
ments. This can be performed through many methods. One 
possible approach to conducting step 112 is to solve the 
known Frenet-Serret equations, as is described in detail in 
Moore ’599 and summarized below. 

Functional Operation 

A key benefit of the present method 100 is that, in the case 
of symmetric core arrangements, it is not necessary to calcu- 
late bend direction in order to calculate curvature. The present 
method 100 is less sensitive to noise at lower strain levels by 
avoiding division by small numbers, whereas in prior art 
techniques, estimation of bend angle can be difficult when the 
strain measurement of a core used in the calculation is small. 

The method 100ofFIG.2canbe executed by the controller 
26 of FIG. 1 using any of a wide variety of commercially 
available hardware platforms running mathematical software 
packages or custom-written code. The method 100 is plat- 
form-agnostic and thus does not rely on a given software 
platform. A multi-core fiber optic cable is connected to an 
interrogation system capable of measuring strain using sen- 
sors embedded in each of the fiber optic core, for instance 
FBGs. 

To execute the method 100, the cable 10 of FIG. 1 must be 
secured at some known point. This represents the initial con- 
dition noted above. Uiis secure point maintains a known 
initial tangent vector. Any shape sensing segment of the fiber 
cable may be free to twist as its own elasticity dictates during 
bending, as would be the case if it were held inside an adhe- 
sive-free tubing, e.g., the sleeve 23 noted above with refer- 
ence to FIG. 1A. Peripheral to the method 100 is a suitable 
method of determining the actual curvature vectors of the 
arbitrary initial shape. As an example, if sleeving is woven 
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into the fabric of an inflatable structure for the purpose of 
containing a multi-core shape sensing fiber, the location of the 
sleeving when the structure is hilly inflated can be indepen- 
dently determined and initial curvature vectors can then be 
calculated for a multi-core fiber that is inserted into the sleev- 5 
ing. 

The accuracy of the curvature and bend estimates con- 
ducted a part of the present method 100 is increased as the 
cores 12 are spaced farther from the center of the shape 
sensing fiber 15 of FIG. 1A. The accuracy is also increased by 10 
the addition of cores 12 to the fiber 15. A four-core shape 
sensing fiber, for example, would tend to be more accurate 
than a three-core fiber, and so on. 

Maintaining the initial condition of the fiber 15 is impor- I5 
tant. Breakdowns of adhesives or other loosening of the initial 
condition mounting arrangement can affect system accuracy. 
The accuracy of the shape measurement will be affected by 
any error in the curvature correction process. If the curvature 
and bend direction of the initial shape contains error, the 20 
accuracy of the subsequently measured shape will be 
affected. 

The prior art of Moore ’599 discloses how to derive a shape 
solution using a multi -core fiber based on bending measure- 
ments made along the length of fiber, which is required to be 25 
free of extrinsic forces. Generally, the Frenet-Serret formulas 
and the corresponding fiber position relationship can be 
solved numerically after explicit functions of curvature and 
torsion are found. Curvature and torsion functions for the 
fiber are found by first calculating discrete sets of curvature 30 
and bend direction at specific points along the fiber, usually at 
sensor locations such as FBGs. In order to deduce curvature 
and bend direction in three dimensions, at least three separate 
core strains must be measured. 

FIGS. 4A and 4B detail the cross-sectional layout of a 35 
typical symmetric tri-core fiber. The locations of the cores 12 
in the cross section of the fiber are defined by a radial offset 
from the center, r, and an angular separation of 2jt/3, with 0j 
representing the offset of core 1 from the 0=0 axis. Not shown 
are the cladding boundary and fiber coating described above 40 
with reference to FIG. 1A. 

The orientation of the y-z axes in the fiber 15 would typi- 
cally align with the orientation of the global y-z axes when the 
fiber 15 is in a straight configuration. The x-axis is pointing 
along the fiber 15 in the direction of increasing fiber length. A 45 
calibration of the orientation of the cores relative to the y-z 
axes is necessary before the calculations of bend direction can 
be carried out. Such a calibration typically consists of bend- 
ing the fiber 15 in the direction of a known plane, measuring 
the induced strain in the bent sensors 24 (see FIG. 1), and 50 
calculating the core orientation which gives proper bend 
direction. The discrete sets of local curvature and bend direc- 
tion can be calculated using known equations, and repre- 
sented as: 


bend angle: ip„=o 

At the “zero” point along the fiber 15, the typical assump- 60 
tion is that the fiber 15 is secured in a known position, usually 
straight, and therefore curvature is assumed to be zero. The 
bend angle at that point is arbitrary and is typically set to zero. 
The “zero” point along the fiber 15 is also usually defined as 
the point at which fiber length is zero: 65 


12 

The discrete local curvature and bend direction sets are con- 
verted into curvature and bend direction functions using a 
curve fitting algorithm: 

curvature: K(.s) 0 . vAn=o v! 

bend direction: ij>(.s) N ,S tl 0 v ] 

The torsion function is found by differentiating the bend 
direction function: 

T(i)=t'(i) 

Initial conditions of the Frenet-Serret frame and fiber position 
are specified at s=0. The initial condition for position at fiber 
length zero is dictated by the actual physical position of the 
fiber in the global reference frame. The trivial choice is: 

(Wo.zo)=(0,0,0). 

The tangent vector initial value. T 0 , is also dictated by the 
actual physical position of the fiber at (x 0 , y 0 , z 0 ), in that T 0 
must be “pointing” in the direction of increasing fiber length 
at s=0 in the global reference frame. A trivial choice is along 
a single axis; such is the case for the data presented in this 
paper, where T 0 is configured to run along the x-axis: 

T 0 =x 

The specification of the normal vector initial value, N 0 , is 
made in the global reference frame and acts to relate the fiber 
local reference frame to global reference frame at s=0. 
Because the normal vector points in the direction of bending, 
specifying N 0 is easily done using: 

Nq = co s (i|i (0) )jp+s in(ij) (0) )£. 

The Binormal vector initial value, B 0 , is found using the 
Frenet-Serret frame definition that Binomial, normal, and 
tangent give an orthonormal frame: 

B 0 =TqxNq 

With initial conditions specified, the Frenet-Serret equations 
can be solved using numerical methods to give location of the 
fiber vs. fiber length, as well as the tangent, normal and 
binormal vectors along the fiber 15. The explicit system of 
equations in detail is: 

T' x (s)=k(s)N x (s) r y (s)=K(s)N y (s) T' x (s)=k(s)N x (s) 

N' x (s)=K(s)TJs)+T(s)B x (s);N' y (s)=-K(s)r y (s)+x(s) 

B y ( S );N' x (s)=-K(s)T x ( S )+T(s)B z (s) 

B' x (s)=-t(s)N x (s ) B' y (s)=-x(s)N y (s) B' z (s)=-x(s)N z (s) 

xXs)=T x (.s); y’(s)=T y (s); z'(s)=TJs) 

Initial conditions are: 

*(0)X0)40),r x (0),fy(0),r z (0)^(0)^(0)^(0)A 

(.0)J3 y (0)J3JO) 

The solution gives the position of the fiber 15 in three-dimen- 
sions as a function of the fiber length: 

Extrinsic Forces 

Conventional methods exist for using multi-core fiber 
bending measurements in the solution of Frenet-Serret equa- 
tions in the absence of external or extrinsic forces. Such 
forces can act on a fiber in such a way as to give measurement 
results that differ from the actual shape of the fiber. One such 
extrinsic force is axial strain. Under axial strain, sheathing of 
a fiber or an adhesive securing the fiber to a surface transfers 
loads into the fiber in the axial direction. While certain intrin- 
sic behavior of the fiber is to be expected for any given fiber 
shape, any extrinsic behavior not directly related to the fiber 


^o=to=-So=0 
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shape can lead to inaccurate fiber shape measurements. The 
example method 200 of FIG. 6, which is described below in 
conjunction with FIG. 5, may be used to incorporate mea- 
sured extrinsic effects into the calculation of shape and/or end 
position of a multi-core fiber, e.g., of the type used in the cable 5 
10 of FIG. 1. 

While a frictionless-sheathed fiber will typically be free of 
externally-induced loads, axial loading of the fiber can be 
externally-induced through attachment mechanisms, e.g., 
attachment to an undersea tow cable, tight weaving in an 10 
inflatable structure, attachment effects from tethered satel- 
lites, and extreme handling by a surgeon in the case of mini- 
mally invasive surgical catheters. In tight -bending situations, 
sleeving that is otherwise frictionless can be elliptical in L - 
cross-section to the point of tightening onto the fiber and 
transferring external loads into the fiber. Thus, the present 
method 200 accounts for extrinsic axial strain behavior of the 
fiber, unlike prior art approaches that assume such axial strain 
is zero. 20 

Referring to FIG. 5, each core 12 is located a distance r, 
from the center of the fiber and has an angular offset from the 
y-axis of 0 ( . The y-z axes represent the local natural orienta- 
tion of a Frenet-Serret frame. The angle €>,, is the angular 
offset of the neutral bend axis from the y-axis. A fiber sub- 25 
jected to bending with no extrinsic axial forces will have zero 
strain along the bend axis of the cross section and the fiber 
length is assumed to be constant no matter the shape of the 
fiber. Conventional approaches such as that disclosed in 
Moore ’599 can be used to deduce shape in such an instance. 3< 1 

However, if a fiber experiences axial forces, the fiber will 
have non-zero strain along its bend axis. In this case, the fiber 
length actually changes. The fiber will lengthen for positive 
strain and shorten for negative strain. The axial effects the ,,- 
fiber experiences can be location-dependent. That is the fiber 
may be experiencing positive strain in some sections and 
negative strain in other sections. The extrinsic axial effects 
which the fiber experiences are dependent on the application 
and can change from measurement to measurement. Various 40 
measurement techniques may be used to determine the 
applied strain along the neutral bend axis of a fiber, as is 
known by those having ordinary skill in the art. See for 
instance the applied strain value e a in equation ( 0 ) above. 

Incorporating axial strain measurements into the solution 45 
of shape may be accomplished by using the axial strain infor- 
mation in the speed of the curve in the Frenet-Serret equa- 
tions. In Moore ’599, curvature and torsion are typically 
defined with the assumption that the measurements are made 
at some spacing that is constant throughout the fiber, such as 50 
at FBGs that are some repeated distance apart from one to the 
next. The method 200 of FIG. 6 uses the measured axial strain 
of the multi-core fiber 15 of FIG. 1 A in determining the speed 
of the curve in the Frenet-Serret equations. This serves to 
shorten or lengthen the fiber 15 as necessary as the equations 55 
are being solved, thus giving the correct results even though 
the fiber length may change from measurement to measure- 
ment. 

As noted elsewhere above, prior art approaches assume 
that an optical fiber is a linear elastic tube with no extrinsic 60 
forces acting in the axial or angular direction. Shape is 
induced into the multi -core fiber of the prior art through radial 
forces acting through a sleeving with an assumed frictionless 
inner surface. Such a frictionless inner surface is not required 
with the present approach. For instance, friction may be 65 
present that can induce twisting of the fiber. Due to such 
assumptions, prior art approaches assume that the material 
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frame of the fiber follows its natural frame, i.e., the Frenet 
frame, of the curve theory describing its shape in three dimen- 
sions. 

By way of contrast, in the method 200 disclosed below, any 
angular forces acting to twist the fiber of the cable 10 shown 
in FIG. 1 are accounted for by including the deviation of the 
material frame from the natural frame in a torsion function. 
The redefined torsion function then includes two angular 
rates of change one describing the rate of angular change of 
the Normal and Binomial vectors as a function of fiber length, 
and the other describing the rate of change of the orientation 
between the material frame and the natural frame. 

Referring to FIG. 6, the method 200 begins at step 201, 
wherein the equations for curvature K(t), bending if (t), twist- 
ing (j>(t), and strain E(t) are generated by curve fitting the 
corresponding discrete measured sets to give the following 
explicit functions: 


*(4?), Wt 0 ), S'), and E(C S'), 

and to thereby obtain the following equations: 

K(f)=Fit{K_o... w f} 

■q*(r)=Fit{ip /=0 N ,t} 

<KO=Fit {<|> / _ 0 . jv.r} 

£(0=Fit{E_ o ... M f} 

At step 203, the controller 26 of FIG. 1 derives a function 
of torsion vs. fiber length. The function may be represented 
as: 

T<r)=*|>’(r)+iK 0 

where if(t) is the twist rate, e.g., in radians/meter. Those of 
ordinary skill in the art will appreciate that twist measure- 
ments can take alternative forms. Specifically, twist rate or 
absolute twist can be used by the controller 26 to derive the 
function of torsion vs. fiber length. Absolute twist may be 
represented as <1> M J), . . .N (radians). Absolute twist and twist 
rate are related, by the following equations: 

<K=f<l> 

Thus, step 203 may be executed by the controller 26 in a 
particular embodiment by differentiating the binding func- 
tion if (t), differentiating the absolute twist function 0(1), and 
then adding the two differentiated functions, i.e., if'(t) and 
<f>'(t) to derive the function of torsion vs. fiber length. This is 
represented mathematically as 

•t(r)=i|)'(9+ < i > '(2) 

At step 205, the controller 26 sets the necessary initial 
conditions r 0 , T 0 , N 0 , and B 0 =T 0 xN 0 for the known Frenet- 
Serret equations. The general form of the Frenet-Serret equa- 
tions to be solved is: 

r(t)=sm+m)KW(t) 

Bxt)=sm+E(mmt) 

r(t)=r 0 +{sm+E(t))T(t)dt 

Initial conditions of the Frenet-Serret formulas and fiber posi- 
tion are specified at t=0. The initial condition for position at 
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fiber length zero is dictated by the actual physical position of 
the fiber in the global reference frame. The trivial choice is: 

r 0 = (xo,yo,Zo)=(0,0fi) 

The tangent vector initial value, T 0 , is also dictated by the 
actual physical position of the fiber at (x 0 , y 0 , Zg), in that T 0 
must be “pointing” in the direction of increasing fiber length 
at t=0 in the global reference frame. A trivial choice is along 
a single axis, where T 0 is configured to run along the x-axis: 


The specification of the normal vector initial value, N 0 , is 
made in the global reference frame and acts to relate the fiber 
local reference frame to global reference frame at t=0. 
Because the normal vector points in the direction of bending, 1 5 
specifying N 0 is easily done using: 

7V 0 =cos(i|> (0))j>+sin(i|> (0))f 

The Binormal vector initial value, B 0 , is found using the 
Frenet-Serret frame relationship that the Binomial, normal, 
and tangent vectors form an orthonormal frame: 

B 0 =T 0 xN 0 

With initial conditions specified, the Frenet-Serret equations 
can be solved using numerical methods to give location of the 
fiber vs. fiber length, as well as the tangent, normal and 25 
Binormal vectors along the fiber. 

At step 207, with s(t) defining the sensor locations or 
measurement locations as a function t, the following system 
may be solved numerically via the controller 26 of FIG. 1: 

30 

rjfirmv+EmwM 

N' x (t)=s'(t)(l+E(t))-K(t)T x (t)+x(t)B x (t)]. 

B ’ x (t)=-s'(t)(l+E(t)yz(t)N x (t), 

35 

r/t^ma+EmwAt), 

^(r)=i'(f)(l+£(0)-ic(r)^(r)+T(r)^(f)), 

b ’ J ,(.t)=-s'm+E(.mt)Ny(.t), 40 

r z (t)=s'm+Emw x (t), 

B'M^rm+EmWM 45 

x‘(t)=s~m+E(t))T x (t), 

y(t)=sm + E(t))T y (t), 

50 

zxt)=mi+m)rM 

The solution to the equations of step 207 gives the position 
{x(t), y(t), z(t)}, the tangent vector {T T (t), T,.(t), T_(t)}, the 
normal vector {N T (t), N (t), hL(t)}, and the bi-nomial vector 
{B T (t), B v (t), B r (t)}. The variable t is a parametric variable 55 
corresponding to a measurement number along the fiber 15 of 
FIG. 1A. A suitable control action can then be taken using the 
determined shape and position, e.g., with respect to the device 
29 ofFIG. 1. 

Hie method 100 and 200 described above with reference to 60 
FIGS. 2 and 6, respectively, thus improve upon prior art 
approaches which assume that a fiber is a linear elastic tube 
with no extrinsic forces in the axial or angular direction. 
Shape is induced into the multi-core fiber of the prior art 
through radial forces acting through a sleeving with an 65 
assumed frictionless inner surface. Due to those assumptions, 
prior art approaches assume the material frame of the fiber 


follows the natural frame, i.e., the Frenet frame, of the curve 
theory describing its shape in three dimensions. 

In the approach disclosed herein, any angular forces acting 
to twist the fiber are accounted for by including the deviation 
of the material frame from the natural frame in the torsion 
function. The redefined torsion function then includes two 
angular rates of change: one describes the rate of angular 
change of the Normal and Binomial vector as a function of 
fiber length and the other describes the rate of change of the 
orientation between the material frame and the natural frame. 

While the best modes for carrying out the invention have 
been described in detail, those familiar with the art to which 
this invention relates will recognize various alternative 
designs and embodiments for practicing the invention within 
the scope of the appended claims. 

The invention claimed is: 

1. A method of determining the shape of a multi-core 
optical fiber, the method comprising: 

initializing a controller, including: 

positioning the multi-core fiber in an arbitrary initial 
shape; 

determining an initial curvature vector of the fiber that 
describes the arbitrary initial shape; and 
initializing strain measurements using the controller; 
and 

determining the shape of the fiber, including: 

positioning the fiber in a desired shape to be measured; 
measuring strain over a length of the fiber using a plu- 
rality of strain sensors while the cable is held in the 
shape to be measured, wherein the strain sensors are 
distributed along a length of each of the cores; 
defining a three-coordinate p-vector for each of the cores 
as a function of the distance of each respective core 
from a center point of the fiber and a bending angle of 
each respective core with respect to a local axis of the 
fiber, wherein the coordinates of the p-vector define a 
plane; 

calculating, via the controller, an applied strain value of 
the fiber using the p-vector and the measured strain; 
calculating strain due to bending as a function of the 
measured strain and the applied strain value; 
defining an apparent local curvature vector for each core 
of the fiber as a function of the calculated strain due to 
bending; 

calculating curvature and bend direction of the fiber 
using the apparent local curvature vector and the ini- 
tial curvature vector; and 

determining the shape to be measured via the controller 
using the calculated curvature and bend direction. 

2. The method of claim 1, further comprising: 

executing a control action with respect to the fiber using the 

determined shape. 

3. The method of claim 1, wherein positioning the fiber in 
an arbitrary initial shape includes positioning the fiber in a 
non-straight initial shape. 

4 . The method of claim 1 , wherein determining the shape of 
the fiber includes accounting for at least one extrinsic force 
acting on the fiber. 

5. The method of claim 1, wherein measuring strain 
includes measuring a corresponding Bragg wavelength of 
each of a plurality of Fiber Bragg Gratings at different axial 
locations along a length of each of the cores. 

6 . The method of claim 1 , wherein determining the shape of 
the fiber via the controller includes using the controller to 
solve Frenet-Serret equations. 
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7. The method of claim 1, further comprising: 

recording discrete measured sets for each of: curvature, 

bending, twisting, and strain of the fiber; 

performing a curve-fitting operation on each measured set 
to generate a curvature function, a bending function, a 5 
twisting function, and a strain function. 

8. The method of claim 7, further comprising: 

using the controller to calculate a function of torsion vs. 
fiber length, including differentiating the bending func- 
tion and then adding the twisting function. 1 

9. The method of claim 8, wherein determining the shape of 
the fiber via the controller includes solving a set of Frenet- 
Serret equations using each of the curve-fitted sets to thereby 
derive the position of the fiber in a three-dimensional space. I5 

10. The method of claim 9, wherein solving the set of 
Frenet-Serret equations includes determining the position, 
tangent vector, normal vector, and bi-normal vector of the 
cable with respect to the three-dimensional space. 

11. The method of claim 1, further comprising: 20 

attaching a cable containing the fiber to a device; 

wherein determining the shape of the fiber includes calcu- 
lating the position of the device in three-dimensional 
space. 

12. A system for determining a three-dimensional (3D) 25 
shape of a multi-core optical fiber, the system comprising: 

a fiber optic cable having an optical fiber with a plurality of 
cores, and further having a plurality of strain sensors, 
wherein the strain sensors are configured to measure 
strain at different axial locations of each of the cores; and 3< 1 
a controller in communication with the plurality of strain 
sensors, wherein the controller is configured to: 
calculate, using the measured strain, an applied strain 
value as an axial strain component experienced by 35 
all of the cores due to at least temperature and an 
axial load acting on the cable; 
calculate strain due to bending of the fiber as a func- 
tion of the measured strain and the applied strain 
value; 40 

define an apparent local curvature vector for each core 
as a function of the calculated strain due to bend- 
ing; 

calculate curvature and bend direction of the fiber 
using the apparent local curvature vector and the 45 
initial curvature vector; 

generate a curvature function, a bending function, a 
twisting function, and a strain function by perform- 
ing a curve-fitting operation; 


calculate a function of torsion vs. fiber length, at least 
in part by differentiating the bending function and 
adding the twisting function to the differentiated 
bending function; and 

determining the shape of the cable via the controller in 
part by solving a set of Frenet-Serret equations 
using each of the curvature function, the bending 
function, the twisting function, and the strain func- 
tion to thereby derive the shape of the cable in a 
three-dimensional space. 

13. The system of claim 12, wherein the controller is con- 
figured to calculate the applied strain value using a p-vector 
and the measured strain values for each core. 

14. The system of claim 12, further comprising: 
a device connected to one end of the cable; 

wherein the controller is configured to determine the posi- 
tion in three-dimensional space of the device using the 
determined shape of the cable. 

15. The system of claim 12, wherein the arbitrary initial 
shape is a non-straight shape. 

16. The system of claim 12, wherein the cable includes a 
sleeve which circumscribes the fiber, and wherein the sleeve 
is at least partially bound to the fiber. 

1 7 . A method of determining the shape of a multi -core fiber 
optic cable in the presence of extrinsic twisting and stretching 
forces, the method comprising: 

recording discrete measured sets for each of: curvature, 
bending, twisting, and strain of the cable; 
performing a curve-fitting operation on each measured set 
to generate a curvature function, a bending function, a 
twisting function, and a strain function; 
using a controller to calculate a function of torsion vs. fiber 
length, including differentiating the bending function 
and then adding the twist function; and 
determining the shape of the fiber via the controller in part 
by solving a set of Frenet-Serret equations using each of 
the curvature function, a bending function, a twisting 
function, and a strain function to thereby derive the 
shape of the fiber in a three-dimensional space. 

18. The method of claim 17, wherein solving the set of 
Frenet-Serret equations includes determining the position, 
tangent vector, normal vector, and bi-normal vector of the 
cable with respect to the three-dimensional space. 

19. The method of claim 17, wherein the fiber is bound to 
a sleeve of the cable such that the extrinsic forces include a 
twisting force imparted to the fiber by the sleeve. 

20. The method of claim 17, further comprising: 
winding the cable onto a spool such that the initial shape of 

the cable is helical. 



